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The Luria-Delbriick experiment is a cornerstone of evolutionary theory, demonstrating the ran¬ 
domness of mutations before selection. The distribution of the number of mutants in this experiment 
has been the subject of intense investigation during the last 70 years. Despite this considerable ef¬ 
fort, most of the results have been obtained under the assumption of constant growth rate, which is 
far from the experimental condition. We derive here the properties of this distribution for arbitrary 
growth function, for both the deterministic and stochastic growth of the mutants. The derivation 
we propose uses the number of wild type bacteria as the independent variable instead of time. The 
derivation is surprisingly simple and versatile, allowing many generalizations to be taken easily into 
account. 


I. INTRODUCTION. 

The Darwin-Wallace theory of evolution rests upon 
mutation of living organisms and their selection. In their 
seminal article [1], Luria and Delbriick (LD) described an 
experiment demonstrating the randomness of mutational 
events before the selection process. The experiment con¬ 
sists of growing C cultures of bacteria in parallel in iden¬ 
tical environments, beginning with a small number Nq 
(typically 10 3 ) in each batch. After a sufficient growth 
period, the cultures saturate and the number of wild type 
bacteria reaches N (typically 10 9 — 10 10 ). Each culture is 
then tested against an antibacterial agent, a phage virus 
in the LD case, and the number of surviving bacteria aris¬ 
ing from mutation in the cultures is counted by a plating 
method. If C is large, the probability P(m) of having m 
mutants can be experimentally determined; in practice, 
C cannot be large and therefore only statistical quanti¬ 
ties such as the mean and the variance of the number of 
mutants can be estimated. 

The LD experiment has spurred a large interest and 
many authors have developed increasingly refined mod¬ 
els to estimate statistical properties of the random vari¬ 
able X of the number of mutants, such as its cumu- 
lant/probability/moment generating functions (cgf, pgf, 
mgf), from which the mutation rate or probability can 
be estimated. The pioneering authors were Lea and 
CoulsonJU, Armitage Q, Bartlett0|, Crump and Hoel[5], 
Mandelbrotj(|, Sarkar, Ma and SandriQ who set the LD 
distribution on a solid mathematical ground, generalized 
the model to take into account stochastic growth of the 
mutant and the wild type, and developed algorithms to 
estimate the mutation rate. The works of these and 
other authors have been reviewed in an elegant article 
by Zhengjgj which also contains original results and cor¬ 
rections of some of the errors contained in the previous 
works. These investigations have been extended during 
the last 15 years by authors such as Angerer[9|, Dewanji 
et al|T(|, Ycart jTl] , Kessler and Levinefl2i.fl3]. A descrip¬ 
tion of some of these more recent works will be given in 
the following sections. 

The fundamental LD experiment is now currently used 


to estimate mutation rates in various setups such as an¬ 
tibiotic resistance or experimental investigation of the 
evolutionary mocess [14l4l6| . 

However, nearly all of the existing computations have 
focused on the exponential growth (either deterministic 
or stochastic) of the wild type and mutant bacteria, al¬ 
though Dewanji et al.(Tol| have extended these results to 
Gompertz growth. The reason behind this choice is that 
in these formulations, an explicit expression for the num¬ 
ber of wild type (WT) cells as a function of time is needed 
in order to compute the statistical properties of the num¬ 
ber of mutants. 

The assumption of constant growth rate is however too 
restrictive. Experimentally, the growth is never exponen¬ 
tial but follows a Monod curve[l7| : the growth rate is 
not constant, but begins with a value close to zero (called 
the lag phase), increases gradually to a maximum value 
and then decreases as the number of bacteria increases, 
to finally reach zero when the culture is saturated. Var¬ 
ious functions (logistic, Gompertz, Richards, Stannard, 
...) are used in the literature to model the growth curve 
and their relevance has been studied in depth by Zwi- 
eteting et al. [Tsl] . 

The real time however is not the relevant independent 
variable in terms of which the system may be described. 
The WT population grows from an initial number of cells 
No to reach a final value N. Each time a WT cell divides, 
there is a small probability v that a mutant having the 
desired trait (phage or antibiotic resistance) appears. It 
does not matter how much time the system spends be¬ 
tween WT population size n and n +1, but only the fact 
that once a division has taken place, a mutation may ap¬ 
pear. For the mutants growing in the same environment 
as the WT, their growth curve will be similar (but not 
necessarily equal) to that of the WT. The only quanti¬ 
ties that are indeed measured in an experiment are the 
number of mutants to, the initial WT population Nq and 
the final population N + to. Even though the growth 
curve can theoretically be measured, its determination is 
cumbersome and, as we will see, not relevant. 

In this paper, we shall use the WT population size n 
as the independent variable. It appears that this formu- 
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lation of the LD distribution is surprisingly simple and 
applies to any growth curve for the bacteria, including 
of course exponential growth. As in the case of classical 
derivation of the LD distribution, this derivation is valid 
in the limit of large final population size, which is typi¬ 
cal of experiments involving bacteria. The formulation is 
versatile and can take into account many generalizations, 
some of which are considered in this article. A similar 
approach has been taken by Kessler and Levine [13, H3| 
when solving directly for the Master equation governing 
the dynamics of the mutant population. 

This article is organized as follows: in the next section, 
we present the basic concept for the simple case of deter¬ 
ministic growth of both WT and mutant bacteria, where 
only mutations apparitions are random. We present some 
generalizations, such as different growth rates for WT 
and mutants and non-constant mutation probability. In 
the following section, we generalize the model to the case 
of stochastic growth of the mutant, where we consider (i) 
a linear birth process for the mutant and (ii) a random 
relative growth rate for the mutant. We stress that in 
the following computations, the growth rate is not con¬ 
stant, but can have an arbitrary form. The last section 
is devoted to a discussion of possible extensions of this 
work and conclusions. An appendix, containing straight¬ 
forward mathematical derivations is included in order to 
make this article self sufficient. 


II. DETERMINISTIC MODEL. 

A. Equal growth of WT and mutant. 

Consider a culture of WT bacteria growing from size 
Nq to size N. The growth curve can be as general as pos¬ 
sible assuming that no death event takes place. Let X. n be 
the random variable describing the occurrence of a mu¬ 
tant when the WT population increases from nton+1. 
Throughout this paper, we use the term mutant to des¬ 
ignate an individual which acquire a trait (e.g. phage or 
antibody resistance) which will be tested once the growth 
is stopped. 

Denoting the mutation probability by v, 

Pr(X n = 0) = 1 - v ; Pr(X n = 1) = u (1) 

This may seem to be an approximate description, be¬ 
cause if during a cell division, a mutation has occurred, 
obviously the number of WT cells has not increased 
from nton+1. A more precise formulation would be 
Pr(A' n = k) = (1 — v)v k , i.e., the number of mutants 
when the WT population increases by one unit is geo¬ 
metrically distributed. However, as v -C 1 (v is usually 
of the order of 10 -8 ), we will use the relation m to de¬ 
scribe the random variable X n . The generalization to 
the geometric distribution of X n is straightforward (see 
appendix QJ . Note that most formulations of LD distri¬ 
bution use the above approximation. 



Figure 1: The number of mutants X at WT population size N 
is the sum of the contributions of mutants lineages appearing 
at WT population size rn. The size of the lineage of mutant i 
appearing at WT population size rn is rrii{n ) = n/rii for n > 
rii. The i —th mutant, appearing at size 1 at WT population 
rii will be present at size N/m in the final population. 

We assume in this subsection that both mutant and the 
WT population follow a deterministic, equal growth. We 
do not assume the growth rate to be a constant. As the 
mutants are similar to the WTs, a mutant appearing in 
one copy at WT population size n will contribute N/n to 
the number of mutants when the WT population reaches 
size N (figure [1]). In other words, the proportion of the 
number of this mutant to the number of WT population 
will remain constant. 

Let Y* = ( N/n)X n be the contribution of this mu¬ 
tant to the final number of mutants X, when the WT 
population reaches size N. Then 

N 

X=J2 Y n 

n=No 

As the mutant occurrences are independent random vari¬ 
ables, the moment (mgf) and cumulant generating func¬ 
tions (cgf) are 

N 

m = <'■*>= n 

n=N 0 

N 

V’(s) = log«e sA ))= J2 log ((e 8 ^)) 

n=No 

On the other hand, by its very definition, 

(e aY ?) = \-v + ue aN / n 

which gives the cgf as 
N 

i/j(s) = log (l - v + ^e sJV/n ) 

n=N 0 

= N f log (l - i/ + ^e s/x ) dx (2) 
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where in the second expression, we have used the con¬ 
tinuous approximation for the sum, x = n/N and Xq = 
No/N. The relative error in using the continuous approx¬ 
imation is at most O (1/Nq log(a;o)). 

Note that this derivation is analogous to the filtered 
Poisson process derivation used when the problem is for¬ 
mulated in real time (Q,eq. 14). However, because the 
problem is formulated in terms of WT population size, 
the propagator is simply a straight line regardless of the 
growth rate function (figure [T|). 

Expanding the expression © to the first order in v 
and restricting the domain of definition to s < — Xq log u, 

= -6(j) + 0 f e s / x dx + 0(v 2 ) (3) 

J Xq 

Where 6 = Nv and (j) = 1 — xq . For the particular case 
of exponential growth, the expression © for the cgf has 
been obtained by Crump and Hoel (S] and in closed form 
by Zheng (||,eq.l4). Note that for Zheng, N = exp{fh), 
i.e. the initial number of WT bacteria is 1. 

The first two cumulant coefficients k p = ^)^(O) are 
then 



Figure 2: Contribution of mutants with different growth rate. 
The size of the lineage of mutant i appearing at WT popula¬ 
tion size Hi is mi(n) = ( n/rii) c for n > m. 

the population size is n will contribute ( N/n) c to the final 
number of mutants (figure © . The computation of the 
preceding section can be repeated and leads to 


«i = — fllogxo 

6(j) 

k 2 = — 
x 0 


(4) 

(5) 


These expressions for the average (k± ) and the variance 
(«2 ) have been obtained originally by LD for the ex¬ 
ponential growth of bacteria [l|. As we see here, the 
hypothesis of exponential growth is superfluous and the 
expressions m are valid for arbitrary growth curves. 
To the leading order in v and Xq , the general expression 
for cumulant coefficients is 


Kp +1 = p P > 0 

P x 0 


( 6 ) 


r 1 

V’(s) =N log (l — v + ve s ^ x \ dx 

' Xo ' ' 

Keeping only the leading term in v and xq , for c ^ 1/ p , 
the p —th cumulant coefficient is given by 

The case c = 1/p can be recovered from the above for¬ 
mula by taking the limiting value for c —> 1/p and reads 

k p = -9 logzo 


These expressions are known for the special case of ex¬ 
ponential growth (js||,eq.9 for equal growth rate). 


For the particular case of exponential growth (a(n, t) = a 
), this is the expression given by Zheng ([s||, eq.9). 


B. Different growth of WT and mutant. C . Variable mutation probability. 


Let us now consider the case where WT (n) and mu¬ 
tants (m) (once they have appeared) have similar but 
different growth rates: 

dn . . dm . . 

—— = a(n,t)n ; —— = ca(n,t)m 
dt y ’ ' dt K ' 

where c is a constant. We do not specify any particular 
form for the growth rate, but we suppose that the mu¬ 
tant follows the same law as the WT, within a constant 
multiplicative factor. This is the case for example where 
the resources are depleted by the growth of the bacteria, 
and the mutant is inferior to the WT for its duplication. 

Time can be eliminated between the above equations: 
dm/dn = c(m/n). A mutant appearing at one copy when 


In most models the mutation probability v is consid¬ 
ered to be a constant and independent of time, i.e., pop¬ 
ulation size. This is a sound hypothesis when the muta¬ 
tion involves only point-mutations on the chromosome. 
However, traits such as virus or antibiotic resistance may 
involve many point mutations before the trait is func¬ 
tional. Bacteria at the end of the growth process, having 
achieved more divisions, may be more prone to mutate 
to the given trait than bacteria at the early stage of the 
growth. A crude approximation of the above phenomena 
will be a mutation rate that depends on the population 
size v = v(n). The formulation for the number of mu¬ 
tants of the previous section does not suppose a constant 
rate of mutation and the relation© for remains valid. 
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Figure 3: The stochastic propagator K„ of the mutant ap¬ 
pearing at WT population size n. 


The first two cumulants are given by 


tion size N is 

N 

*K*)= E log(l-^ + ^(e sif ")). (8) 

n=No 

Knowing the statistical properties of the propagator gives 
access directly to the statistical properties of the total 
number of mutants. Before applying this concept to spe¬ 
cific cases, let us compute the first two cumulant coeffi¬ 
cients: 

N 

*1 = * E ( K n ) (9) 

n=N 0 

N N 

« = -£ ( W ) 2 )-- 2 £ «> 2 ( 10 ) 

n=No n=No 

The mean k± is what we already had in the deterministic 
case, where (K„) = N/n. Let us express the second 
moment of K x as a function of its mean and variance 
V N 

v n 


K1 


H2 


N 

N 


1 


dx 


f 

Jx 0 


1 v{x) — u 2 (x) 


dx 


III. STOCHASTIC GROWTH OF MUTANTS. 
A. General discussion. 

Until now, we have considered the deterministic prop¬ 
agation of the mutant from its appearance at population 
size n to its final value at population size N. We will 
denote the propagator as K x . The contribution of the 
mutant appearing at WT population size n to the final 
population N was expressed as 

V N_ Y K n 
1 n — A n n n 

For the deterministic case, = N/n and the mgf of 
was simply 


(e sY ") = l-v + ve sN / n 

We will now consider for the mutant a stochastic prop¬ 
agator K// (figure [3]). Because X n takes only the values 
0 or 1, 


(e sY * ) = (e sX " K * ) = l-v + v (e sK " ) (7) 

(see appendix [2]) . Therefore, all the discussion of the pre¬ 
ceding section naturally generalizes to stochastic propa¬ 
gation and the cgf for the number of mutants at popula- 


((K x ) 2 )=V n N + (K X ) 2 

Then 

N N 

k 2 = v(l - v) (K x ) 2 + V V x 

u=Nq u—Nq 


The first term on the RHS of the above relation is what 
we already had in the case of deterministic growth. The 
second term is the contribution of the stochasticity of the 
propagator to the variance of the number of mutants at 
population size N. 


B. Linear birth process. 


Consider the case where the growth of the WT is de¬ 
terministic and continuous 


dn 

dt 


= a(n, t)n 


while the mutant, once it has appeared, follows a stochas¬ 
tic growth with transition probability density 

W{m —>• m + 1) = a(n, t)m 


where W(m —> m + 1 )dt is the probability that this lin¬ 
eage of the mutant has increased its size by one unit in the 
interval [t, t + dt}. This model was first introduced by Lea 
and Coulson(2] for exponential growth case a(n, t) = a. 

A note of caution should be made here. Although 
widely used, the linear birth model may not be very re¬ 
alistic, as bacterial division times are not exponentially 
distributed. Indeed, after a division, a bacterium needs 
to elongate again to its original size before being able to 
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divide again, so that the next division time cannot be 
smaller than a finite time r. In fact, the distribution of 
division times around the time r is fairly narrow and the 
division process is much less random than a linear birth 
process. The phenomenon has been experimentally in¬ 
vestigated by a microfluidic device by Wang et al [19| ; 
the overestimation of the mutation rate by a linear birth 
model, for the exponential growth case, has been inves¬ 
tigated by Ycart [ljj. 

Let us now come back to the linear birth model. As in 
the previous section, the real time is not the best choice 
of independent variable and we can write the stochastic 
growth of the mutant in terms of WT population size: 
noting W{rn —> m + 1 )dn the probability that a mutant 
has divided when the WT population size £ [n,n + dn], 
we have 


W(m —» to + 1) = W(m —> m + 1)— = — (11) 

an n 

Note that the equation for the mean of this lineage is 

d (to) (to) 

dn n 

which conserves the ratio between the size of this lineage 
and the WT population, in agreement with the deter¬ 
ministic case investigated above (subsection III Al) . 

The master equation governing the growth of the mu¬ 
tant is 


growth case can be recovered from the above expression 
and is equal to the expressions given by Lea and Coulson 
and Zheng (Q, eq.52-53). 

Other cumulant coefficients can be readily recovered 
by multiple derivation of expression m- Restricting the 
computations to the leading order of v and Xq , the ex¬ 
pression for the cumulant coefficients is (see appendix 

HI): 


k p = 6 


p\ 


(p ~ l)zj 


—T p > 2 

t p— 1 1 — 


(15) 


Note that even in the case of exponential growth, no 
general expression for the cumulant coefficients could 
be obtained by classical methods ([§f). Comparing the 
above expression with the deterministic case where k p = 
Nv/ip— 1):Eq _ 1 , we see that indeed the linear birth pro¬ 
cess induces large amplification of the p —th cumulant 
coefficient by a factor of pi. 

The probabilities. To compute the probabilities, it 
is more advantageous to use the probability generating 
function (pgf) 


G(z) = (z x ) = e^ z) 

where ip( z ) is defined from (fl4l) by setting 2 = e s , i.e .: 

^ = N I„ >»« (' - “ + , (3 ,-?) + i ) * < 16 > 


dP(m; n) 
dn 


= W(m— 1—> m)P(m — T,n)— 

W (to —y to -(- 1).P(to; n) (12) 


and the mgf of the propagator K x is given by (see ap¬ 
pendix HI) : 



ne s 

(■n — N)e s + N 


(13) 


Approximating the above expression to the first order in 
v, we have 

4>{z) = 0{- - 1) log (1 - <j>z) + 0{v 2 ) (17) 

where 9 = Nv and <j> = 1 — Xo- We therefore obtain a 
simple expression for the probability generating function 

G(z) = (1 - </>z) 0(1/z - 1} (18) 


Using now the expression © and the continuous variable 
x = n/N, we obtain the cumulant generating function of 
the number of mutants 

<«»)=jv £ log (i - -+ e , ( /_ e ; T )+1 ) (w) 

Note that the above integral can be expressed in an an¬ 
alytical, albeit cumbersome, form. The expressions for 
the two first cumulant coefficients are 

Ki = —9logxo 

— — (2 - v) + 9 log xq 
x 0 

96 

~ 2— +9\ogx 0 

X 0 

Where as before, 9 = Nv and <j> = 1 — Xo- The variance 
of the number of the mutants is now approximately twice 
what we had for the deterministic case. The exponential 


For the case of exponential growth, the above expres¬ 
sion (without the <f> factor) was first discovered by Lea 
and CoulsonQ; omitting the <j> factor however results 
in divergent moments. The correct expression can be 
seen in the Zheng review (Q, eq 65). We stress again 
that relation (IT51) is very general and does not depend on 
the assumption of exponential growth. Note that in the 
case of exponential growth, 9 ~ exp (/?£) and all the mo¬ 
ments diverge as t —> oo. This divergence, discussed by 
Bartlett and later by Zheng [8], cannot be cured within 
the framework of the exponential growth model. No such 
divergence exists in the present formulation, as the WT 
population size, following any realistic growth curve, will 
remain finite. 

In order to evaluate the probabilities, we have first to 
compute Expanding expression HID in powers 

of z we have 


1 p{z) 


-9^ + 9j2 

n= 1 



A n +! 


Z 


n 


n 


n + 1 
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Keeping only the first leading terms in 9 and xq we have 


ip(z = 0) 

= 0 ) 
p\ 


—9(f> 

9 9px% 

p(p+ 1) P + 1 


(19) 

( 20 ) 


where the second term in Xq can be neglected for p -C 
l/a;o and 9f> ss 9. We can now use the Faa Di Bruno 
Formula [20j to compute the k —th derivative of G(z) at 
z = 0 and obtain the probabilities : 


P{k) = 


= e 


8 



mj 


where the sum is taken over all k —tuples {mj} such that 
jmj = k. The above formula can be easily pro¬ 
grammed to compute numerically the probabilities. The 
only linear term in 9 in the above formula is for m*, = 1, 
rn^ot = 0. Therefore, for 9 <C 1, the probabilities take 
the simple form of 


m 

m 


„-0 
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k(k + 1) 


k > 0 


The explicit expression for the probabilities in the linear 
birth model has been intensely investigated by many au¬ 
thors, and reviewed by Zheng ([§}, 5.3) and the above ex¬ 
pressions are known for the constant linear birth model. 
An explicit expression in terms of Landau distribution 
has been found recently by Kessler and Levine flij . Il3| 
for both the deterministic and stochastic growth model. 

Experimentally, obtaining the probabilities implies a 
very large number of parallel cultures and the above ex¬ 
pressions may not be of great practical use. 

Different growth rate. The discussion of the preced¬ 
ing subsection can be extended to take into account a 
different relative growth rate for the mutant compared 
to the WT: 



Figure 4: Deterministic growth, random relative growth rate. 
The size of the lineage of mutant i appearing at WT popula¬ 
tion size Hi is rrii(n ) = ( n/ni) c for n > m. This time however, 
c is a random variable. 


C. random relative growth rate. 

The traditional formulation of LD distribution assumes 
that all mutants are similar in their growth function. As 
many different mutations can bring a bacteria to the same 
phage resistance, this assumption may seem too restric¬ 
tive and can be easily relaxed. For example, a compre¬ 
hensive study of this phenomenon has recently been pub¬ 
lished mi where the growth rate of all mutants in the 
gene TEM-1, conferring resistance to the antibiotic ce¬ 
fotaxime, where measured and shown to be variable in 
some conditions. 

Consider the case where the growth of both mutants 
and WT are deterministic as in subsection III Bl 

dn . dm . . 

—— = a(n,t)n ; —— = ca[n,t)m 
dt v ' dt K ' 

but the relative growth rate c is now a random variable 
(figure |U): when it appears, a mutant picks a relative 
growth rate c from a given distribution, which is trans¬ 
mitted to its progeny. 

Following the discussion of subsection IIII Al the prop¬ 
agator this time is 


dn . 

= a(n,t)n 
dt v ’ ' 

W(m—>7n+ 1) = ca(n,t)m 

Repeating the discussion of the preceding sections, the 
cumulant generating function in this case is (see appendix 

[1 

m = N £ log (i - - + - t) < 21 > 

from which the cumulant coefficients can be deduced as 
before. 


K% = (. N/nf 

where c is now a random variable. Let us denote p(s) its 
moment generating function 

p( S ) = (e sc ) 

Then according to relations (liJlfTUl). the first two cumu- 
lants are now, to the first order in v: 
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For example, if c follows a Normal distribution c = 
Af(n,cr), then p(s) = exp (ps + (1/2)cr 2 s 2 ) and these ex¬ 
pressions can be evaluated in terms of the error function. 
For the specific case of p = 1 and a -C 1, restricting the 
computation to the leading orders of v and Xq , we have 

m = —Ologxo ^1 + ^-(logiro) 2 
H 2 = — (l + 2 ct 2 (l + (1 + loga: 0 ) 2 )) 


IV. DISCUSSION AND CONCLUSION. 


In the previous sections we have given the general so¬ 
lution of the LD distribution through the derivation of 
its cumulant generating function. The key point to this 
derivation is to change the independent variable from 
time to WT population size, which is indeed the relevant 
variable : whatever the growth function, a mutation can 
occur only when a WT cell divides and WT population 
size changes from n to n+ 1. This consideration consider¬ 
ably simplifies the solution of the problem and allows us 
to extend the solution to arbitrary growth curves for the 
WT population. The mathematical formulation applies 
equally well to the case of deterministic and stochastic 
growth. 

This mathematical formulation is sufficiently simple to 
allow for many generalizations, some of which have been 
considered in this article: for example variable mutation 
probability (subsection III C'l) or random relative growth 
rate (subsection IIII Cl) have been investigated. 

Other generalizations that we have not developed can 
be considered. For example, for the stochastic growth 
case, we have only considered the linear birth process. 
Other, more realistic cases can be envisaged where the 
distribution of the division times are not exponential, 
along the lines developed by Ycartjll]]. One can also con¬ 
sider the case where both mutants and WT grow stochas¬ 
tically. These generalizations would be straightforward if 
the moment generating function of the propagator K^ 
can be derived in explicit form. Another possible gener¬ 
alization would be to take into account the experimen¬ 
tal uncertainty on the initial and final value of the WT 
population, and its influence on the estimation of mu¬ 
tation probability, as has been considered by Ycart and 
V'eziris|22||. 

To summarize, we have developed in this article a ver¬ 
satile method for investigating the Luria Delbriick distri¬ 
bution with an arbitrary growth function. The method 
uses only very few measurable parameters, namely the 
initial and final number of the WT population. We be¬ 
lieve that the method we propose here can be used as a 
simple basis for further investigations of the LD distribu¬ 
tion. 


Appendix: Mathematical details. 

1. Geometric distribution of mutants appearing at 
WT population size n. 

We have noted v the probability of apparition of a 
mutant when a WT cell divides, and X n the random 
variable tracking the number of mutants appearing when 
WT population changes its size from n to n + 1. The 
probability of producing k mutant during this change is 
geometrically distributed: 

Pr(A„ = k) = (1 - v)v k 

As v <C 1, we have approximated X n by a binary process 
(X n = 0,1) in the article. This constraint can be relaxed. 
Following notations of subsection III Al 

y n\ = \- v 

1 1 - ve sN ! n 

and 

N 

4>(s) = J2 log (( eSY " V )) 

n=No 



where 9 = — AHog(l — v) and <f> = l — Xq. The above 
expression is equal to expression @ to the first order in 
v. 


2. Moment generating function of Y„ in the 
stochastic case. 

Consider the random variable = X n where X n 
is a boolean variable P(X n = 0) = 1 — v and P(X n = 
1) = v ((subsection IIII Al) : For simplicity, we suppose that 
is a positive discrete random variable. Then, 

P(Y n N = 0) = 1 — z' + vP(K™ = 0) 

PiXn =k^0) = VP{K% = k) 

The mgf is then 

OO 

(e SYn ) = J2 P ^ = k ) eSk 

k—0 

= 1 - V + vP{K™ = 0) + 

OO 

5>P(A^ = k)e sk 

k =1 

= l-v + v(e' K Z) 
which is the expression (0. 












3. Moment generating function of the propagator 


A mutant appearing in one copy when the WT popu¬ 
lation size is no will reach size AT" when the WT popu¬ 
lation reaches size n. The forward master equation gov¬ 
erning the probabilities of the propagator, Pr(AT™ 0 = m) 
is given by equation m- The mgf of the propagator, 
</>(s, n) therefore obeys the following equation : 


This equation can be transformed into equation (IA.1I) by 
a simple scaling n —>• n c and the mgf is therefore given 
by 


</>(s,n) 


(n 0 /n) c e s 

((n 0 /n) c - 1) e s + 1 


from which the cumulant generating function of the num¬ 
ber of mutants can be deduced. 


-g + P-Og-O (A.1) 

with the boundary conditions 

4>(s, n = n o) = e s 
4>(s = 0, n) = 1 


4. Cumulant coefficients for the stochastic growth. 

To the first order in u, the cgf for the linear birth model 
(subsection IIII Bl) is given by 


Equation (IA.1I) is a linear first order partial differential 
equation that can be solved by the method of character¬ 
istics: 




(n 0 /n)e s 
(n 0 /n - l)e s + 1 


Changing now the notation to denote n as the WT popu¬ 
lation size of the mutant occurrence, and N the final size 
of the WT population, we obtain the expression (THU) . 

If the relative growth rate of the mutant is not 1 but 
c, where c is an arbitrary constant, the transition proba¬ 
bility for the mutant once it has appeared is 

777 - 

W(m —> m + 1) = c— 
n 


and (f>(s,n ) obeys the following equation: 


1 d( t > , n n 

c n ^ + (1 - e) 5J = ° 


(A.2) 


i^(s) = (e s - 1) log (1 - e s + x 0 e s ) 

Expanding the above function into powers of (1 — e~ s ), 
we have 


V’(s) = (s + logxo)(e s -l) + 

^ i (1 - e ~ s ) n 

/ 77 — 1 r r n ~ ^ 

n=2 U X 0 

Evaluating the p — th derivative at s = 0, we have 


V> p (o) = 


p< 


(P ~ l)a 


— + o(— 2 ) 


which is the expression given in equation (11511 . 
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